1 



Atomic Bose Condensation and the Lattice 

Guy D. Moore a 

a Department of Physics, University of Washington, Seattle, WA 98195, USA 



o 
o 

(N 
Or 



> 

(N 

o 

On 
O 

o 



X 



I show how interaction corrections to the Bose condensation temperature of an atomic gas can be computed 
using a combination of perturbative effective field theory and lattice techniques. 



1. Introduction 

This is a lattice talk, but at first it won't seem 
like it; it will be about atomic systems, and ef- 
fective field theories. But studying the effective 
description of the atomic system will lead to a 
question suited to lattice techniques. I am sum- 
marizing 3 papers ||-|], all with Peter Arnold; 
the middle one contains all the lattice stuff. 

Many systems exhibit Bose condensation; for 
instance, liquid Helium, Cooper pairs in super- 
conductors, and arguably nucleon Cooper pairs in 
large nuclei. These examples condense because of 
interactions; none exhibit Bose condensation as 
originally proposed by Bose and Einstein, where 
condensation occurs exclusively because the tem- 
perature is so low that the only place to store the 
bosons is in the ground state. 

The interest in atomic Bose condensation 
mainly arises because this is what does happen 
here; mutual interactions are not that impor- 
tant, the dominant reason for condensation is 
that there is not enough heat to store the density 
of atoms present in the excited states. In fact, 
interactions must be weak in an atomic Bose con- 
densate; otherwise it would condense into a solid. 

The beauty of small interactions is it means we 
can calculate; the coupling is weak enough that 
perturbatively based "field theorist" techniques, 
like effective field theories and loopwise expan- 
sions, can be used. 

The question I will address in this talk is, how 
do interactions modify the relation between den- 
sity and the temperature of the Bose conden- 
sation transition? For simplicity, and to follow 
the literature, I consider a homogeneous system. 
(This also makes the question well posed.) It 



means that the problem addressed is a little arti- 
ficial and has less contact with experiment than 
we would like; but nevertheless since there is al- 
ready a surprisingly big literature on the problem 
as posed, it is worth addressing. 

2. Effective theories 

The right effective description of a trapped 
atomic gas of bosonic atoms, all in the same spin 
state, is a second quantized Schrodinger equation: 



C = 



d 3 i 



d 3 xd 3 y V(x,y) &V(x)&V(y) . (1) 



The atomic number density is with a sub- 

traction implied so the density vanishes in vac- 
uum. Call this Theory I. Technically, Theory I is 
already an effective theory, containing no infor- 
mation about atomic structure; it is only valid on 
length scales larger than the atomic size. Since 
we are concerned with a dilute system we will 
be interested in scales large compared with the 
scattering length, in which case the nonlocal in- 
teraction term should be replaced with a local 
term plus higher dimension corrections. Further- 
more we are interested in the behavior at finite 
temperature, and will only ask about static prop- 
erties (namely the number density), so we can 
go to the Matsubara formalism, with Euclidean 
periodic time. I work in the grand canonical en- 
semble, so I introduce a chemical potential fi, and 
because the system is homogeneous, V(x) can be 
dropped. I will also set h = 1 from now on. The 

appropriate effective theory is therefore 
rP r 

d xCe , 



Sn = 



dr 



(2) 
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£e = *t(-3 r -^- M W^(*t*) 2 ,(3) 
\ Zm J m 

plus high dimension operators. The number den- 
sity is still ^f^, with its zero temperature value 
subtracted. This effective theory requires regular- 
ization, and I use MS dimensional regularization. 
a sc is the scattering length; we could derive it by 
matching to Theory I, but we rather define it to 
be the coefficient shown. There are also high di- 
mension operators, containing either more fields 
or derivatives; but at the order of interest they 
can all be dropped. I will call this Theory II. 

Theory II contains three interesting length 
scales A: 

1. A ~ a sc ~ 5nm, where it breaks down; 

2. A - A T = y/2n/mT - 100-300nm, the 
thermal scale, where nonzero Matsubara 
frequencies decouple and the behavior be- 
comes effectively 3D; 

3. A ~ A T /a sc , the scale were the behavior be- 
comes strongly coupled. 

The physics of Bose condensation is associated 
with the third scale, where the nonzero Matsub- 
ara modes have decoupled and the physics is ef- 
fectively 3D. Therefore we can write a still sim- 
pler effective theory which contains the physics of 
Bose condensation, which I call Theory III: 

Sm = J d 3 x C m (4) 

i v 
£ 3D = _0(_V 2 + r)0+-(0 2 ) 2 , (5) 

with <j) a two component real field (the two com- 
ponents correspond to the real and imaginary 
parts of ^). Again there are in principle high 
dimension terms, but we won't need them. The 
possibility of such matching was realized by Baym 
et. al. (§. 

The relation between fx, a sc , and (\E'^\[') in The- 
ory II and r, u, and (0 2 ) in Theory III is deter- 
mined by matching between the theories. The 
coupling is weak at A = At, so the matching 
may be done perturbatively. At tree level, u = 
967r 2 a sc /A T , {4> 2 ) = (* t *)/TO, and r = -2mfi. 
The operator <fi 2 mixes with the identity. A one 



loop matching calculation of the mixing gives the 
free theory value for the particle number, 

n bee = C(3/2) At 3 . (6) 

To find the 0(a sc / At) corrections to the num- 
ber density, we need to find this mixing to two 
loops, r to one loop, and the other quantities at 
tree level. We can do even better, and determine 
the 0(a 2 c /A T ) corrections, by going to one higher 
loopwise order in every variable. This is done in 
|| . We cannot go further without including high 
dimension operators, so we will stop at this order. 

Besides the matching calculation, we also need 
the critical values of two quantities in Theory III. 
At linear order in a sc we need the value of {4> 2 ), 
evaluated at the critical value r cr ; t (the critical 
/i). At second order in a sc , we also need to know 
the critical value r cr i t . r cr it is also needed to find 
the total number of atoms in a wide but finite 
trap §. 

3. The lattice 

The values of (<j) 2 } and r cr ;t are nonperturba- 
tive. Theory III is the 3-D 0(2) model, and its in- 
frared behavior falls in the x-y universality class. 
The infrared behavior is strongly coupled. Fur- 
thermore, the quantities of interest, {(f) 2 ) and r cr i t , 
are not universal. They cannot, to our knowledge, 
be reliably determined using known analytic tech- 
niques (perturbation theory, e expansion, high 
temperature expansion, etc.). We turn instead 
to a lattice determination. 

We will refer to the lattice regularization of 
Theory III as Theory IV. (j>, a two component 
real field, is discretized on a 3D cubic lattice of 
spacing ai a tt, with action 

<Siv = afatt X! ^ Latt > ( 7 ) 

£ La tt = ^ L (-^V 2 +r L )0 L + ^(0 2 ) 2 - (8) 

Relating lattice parameters and Theory III pa- 
rameters, including (0 2 ), requires another match- 
ing calculation, which is identical in structure 
to the matching performed between Theories II 
and III; it can be done by lattice perturbation 
theory. In [g] we perform this matching to the 
same, and in some cases higher, loop order, as 
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we used in matching Theories II and III. The 
matching is an expansion in aiatt^, which means 
we want ai at t to be relatively small. However 
we also want to be deep in the critical scal- 
ing regime, requiring a large physical volume. 
To improve the matching to Theory III we use 
an improved lattice Laplacian, containing anti- 
ferromagnetic next-nearest couplings (analogous 
to using the Symanzik rather than Wilson action 
in pure glue QCD). This frustrates the cluster al- 
gorithm, but the multi-grid algorithm || is avail- 
able and proves quite efficient for Monte-Carlo 
use. 

We determine the critical value r cr it by the 
method of Binder |7j] , which accelerates the large 
volume convergence; extrapolation of r cr it to in- 
finite volume is well behaved. Our matching pro- 
cedure leaves linear in aiatt errors in tl, so this 
extrapolation is less trivial: 
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Figure 1. Zero spacing extrapolation for r cr jt- 

The result is that, at renormalization point ~fi = 
u/3, r crit /u 2 = 0.0019201(21) @. 

For ((f> 2 ) the large volume extrapolation is 
tougher, and requires a careful accounting of the 
expected critical behavior. However, the match- 
ing has eliminated more of the lattice spacing cor- 
rection, see Fig. ||. After a double extrapolation, 
we find (0f D )/ti = -0.001198(17) §. 

Putting these results together with the results 
of the matching calculation between Theory II 
and Theory III, we determine the correction to 
T c , at fixed number density n, to second order in 



Oscatt?1 



1/3 



the scattering length. Defining p — 
T c = T Q (1 + c lP + (4 ln(p) + c'i) p 2 ) , (9) 
ci = 1.32(2) , c' 2 = 19.75 , 4 = 75.7(4) , (10) 



where Tq is determined by inverting Eq. (||). 
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Figure 2. Top: volume behavior, fixed ai a tt- 
Bottom: ai a tt behavior, fixed volume. 
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